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Abstract 

A kind of fluid dynamic description for the collective movement of pedestrians 
is developed on the basis of a BOLTZMANN-like gaskinetic model. The differences 
between these pedestrian specific equations and those for ordinary fluids are worked 
out, for example concerning the mechanism of relaxation to equilibrium, the role of 
"pressure" , the special influence of internal friction and the origin of "temperature" . 
Some interesting results are derived that can be compared to real situations, for 
example the development of walking lanes and of pedestrian jams, the propagation 
of waves, and the behavior on a dance floor. Possible applications of the model to 
town- and traffic-planning are outlined. 



1 Introduction 



Former publications on the behavior of pedestrians have been mostly empirical (often in 
the sense of regression analyses) and had the intention to allow planning of efficient traffic 
|T^, ^ Meanwhile there also exist theoretical approaches to pedestrian movement 
0) B 0! B IH) il) © 1131 • However, most theoretical work has been done in the related 
topic of automobile traffic (see, for example, |1|, ^ Hj El' EH)- Especially, there have 
been developed some BoLTZMANN-like (gaskinetic) approaches ^ p^ . 

As far as pedestrian movement is concerned, the author has observed that footprints of 
pedestrian crowds in the snow or quick-motion pictures of pedestrians look similar to 
streamlines of fluids. It is the object of this paper to give a suitable explanation of the 
fluid dynamic properties of pedestrian crowds. Henderson was the first, who compared 
gaskinetic and fluid dynamic models to empirical data of pedestrian crowds |ll2|, [13], ^ 



However, his work started from the conventional theory for ordinary fluids, and assumed 
a conservation of momentum and energy to be fulfilled. In contrast to Henderson's 
approach, this article develops a special theory for pedestrians — without making use of 
the unrealistic conservation assumptions. 

We shall proceed in the following way: First the pedestrians will be distinguished into 
groups of different types /i of motion, normally representing different intended directions 
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of walking. At a time t the pedestrians of each type /i of motion can be characterized 
by several quantities like their place x, their velocity ■u^ and their intended velocity 
(that means the velocity they want to walk with). So, in a given area A can be found 
a density p^^x.v^.iP^^t) of pedestrians having a special type of motion /i and showing 
approximately the quantities x, w^, -u^ at time t. For the densities p^{x, v^, v^, t) equations 
of motion can be set up (section 2). From these equations we shall derive coupled differ- 
ential equations for the spatial density (p^) of pedestrians, their mean velocity {v^) and 
velocity variance {{Sv^^iY) (section 3). The resulting equations show many similarities to 
the equations for ordinary fluids, but they contain some additional terms which take into 
account pedestrian intentions and interactions (sections 3.1, 4.1, 5.1 and 6). In section 4 
we shall treat equilibrium situations and the propagation of density waves. However, in 
lionequilibrium situations the final adaptation time to local equilibrium gives rise to inter- 
nal friction (viscosity) and other additional terms (section 5). Effects of interactions (that 
means of avoiding maneuvers) between pedestrians will be discussed in section 6. These 
effects will lead to some conclusions applicable to town- and traffic-planning (section 7). 

Readers who are not interested in the mathematical aspects may skip the formulas in the 
following. However, the mathematical results are important for analytical, computational 
or empirical evaluations. 



2 Gaskinetic equations 

Pedestrians can be distinguished into different types fi of motion, for example, by their 
different intended directions := "W^/U'W^II of motion (normally 2 opposite directions, at 
crossings 4 directions). More exactly speaking, a pedestrian shall belong to a type /i of 
motion if it wants to walk with an intended velocity 

where 

K ■= K) 

is one of several disjoint and complementary sets. A type n of motion still contains pedes- 
trians with a variety of intended velocities v^^, but the advantage resulting from a suitable 
choice of the sets is to get approximately unimodal densities p^(x, v^, 5^, t) and there- 
fore to obtain appropriate mean value equations (see section |]). p^{x,Vfj,,v'^,t) describes 
the number of pedestrians of type p within an area A = A{x) around place x having 
approximately the intended velocity v^, but approximately the actual velocity v^. More 
exactly, is defined by 

where A^^ is the number of pedestrians of type p which are at time t in a state 

{x',u'^)eU{x)xV{u,) 
out of the neighbourhood W(x) x V(-u^) of x and 

■= (^M' ^ ) • 
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State {x, u^) is an abbreviation for the property 

(x, u^) := (x, v^, t^) , 

that an individual is at place x and wants to walk with the intended velocity t^^, but in 
fact walks with velocity t?^. 

U{x) := {x* eM: \\x* - x\\i < r} (2) 

is a neighbourhood around the place x and belongs to the domain M., which represents 
all accessible (or public) places x. A = A{x) denotes the area of U{x). Similarly, 

V(u,) := {n; = (tr,, : - n^H, < ,5 g U,} 

is a neighbourhood of := (w^, v'p with a volume V = V{Ufj). 

We shall now establish a set of continuity equations, which are similar to the ansatz of 
Alberti and Belli [jl|: 



:= ^^^ + T.S,u + T.C,u + q,. (3) 

These equations can be interpreted as gaskinetic equations (see chapters 2.4, 2.7 of Ref. 
P^ , and §3 of Ref. ^^)- denotes the average mass of pedestrians belonging to type 
p of motion. Apart from special situations it will not depend on p, that means ~ m,j. 
The forces := m^v^ can often be neglected. However, they may be locally varying 
functions, depending on the attraction of the places x. The temporal change iP^ of the 
intended velocity can in principal be a function of place x and time t, but it is normally 
a small quantity ({^ ~ 0). Otherwise, the type p of motion is changed. 

According to (§) the change of the density with time is given by four effects: 

• First by the tendency of the pedestrians to reach their intended velocity "W^ [|T], 0- 
This causes p^ to approach 

plix, v„ vl, t) := 6{v, - ^^)p%x, t^, t) (4) 

(the equilibrium density in the absence of disturbances) with a relaxation time 

7/. 

(see 0). p° is the density of pedestrians with intended velocity -y^ but arbitrary 
actual velocity t/^. 5{.) denotes the DiRAC delta function, which is different from 
zero only where its argument (.) vanishes. 
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Second by the interaction of pedestrians, which can be modelled by a Boltzmann- 
like stosszahlansatz fT^, |18|, |l9l . If we take into account that the interactions are of 
short range (in comparison with r, see (0)), we have: 



Sf^u = III cr^y{%r^, ul; M^, u„; x, t)p^{x, m* , t)py{x, ul, t)d^u^ d'^v*^ d^u^ 

(^fiuiUfj,, Uu] M* , M*; X, t)p^{x, u^, t)pu{x, Uu, t)d^Uu d'^u*^ d^ul . (5) 



This term describes pair interactions between pedestrians of types p and p oc- 
curing with a total rate proportional to the densities p^ and py of both interact- 
ing types of motion. The relative rate for pedestrians of types p and v to change 
their states from {x,u^), {x,Uy) to {x,u*^), {x,ul) due to interactions is given by 
d'fiu{Uf^,Uy;u*^,ul; x,t). Assuming that only the actual velocities -i^i, Vy and not the 
intended velocities v^, are affected by interactions, we obtain 



This results in 

„-?0 +\j2-iO 



S^y{x,v^,t) := S^,y{x,v^,v^,t)d 



with 



(^f^v{v*^^ vt] Vf,, Vy)pf,{x, v*^, t)py{x, vl, t)d Vyd v'^d vl 
(^f^uiVf,, Vy] vl, vl)p^{x, v^, t)py{x, Vy, t) d'^ V y rf^?})^ d'^ vl (6) 
(^lv{v*^^ vt] v^)p^{x, vl, t)py{x, vl, t)d^vl d^vl 
f^/l^(^V' ^'^S ^)Pm(^' ^/^' t)py{x, Vy, t)d^Vy d^ v^ (7) 



:= / (y^y{.,.;.,v)d^v. 



which is used later on. (^) is similar to (|^), and can be interpreted analogously. The 
explicit form of cr*j, will be based on a microscopic model for the interactions and 
is discussed in section 6. 

Third by pedestrians changing their type p of motion to another type ly of motion, 
for example when turning to the right or left at a crossing or when turning back 
(change of the intended direction). This can be modelled by 



C^y{X,Ui„t) = j a1^{Uy;U^;X,t)py{x,Uy,t)d Uy 

cr^'^iu^; Uy] X, t)p^{x, u^, t)d'^Uy 



with a transition rate proportional to the density of the changing type of motion. 

If we assume that for the moment of changing the intended velocity tJjJ the actual 
velocity still remains the same (but of course not thereafter), we have 

o"u'^(mi; ^2; X, t) = al;'^{v^; v^; x, t)6{v2 - Vi) . 
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This results in 
where 



and 



Vi,, t) := j , -i^, -yj, t)o?^'yJ . (10) 



Fourth by the density gain q'^{x,v^,'if^,t) or density loss q~{x,v^,'iP^,t) per time 
unit, caused by pedestrians which enter or leave the system at a marginal place 
X G DM. (for example a house) with the intended velocity := & and the 
actual velocity v^: 



v^, -uj, t) := g+(f , v^, v^, t) - {x, v^, -uj, t) . (11) 



3 Macroscopic equations 



For further discussion we need the notations 



(V'm(^m>'?)> — / V'(^V>'?)- 



gf,{x,v^,t) 

— * 

3lJ.,i 



mf,Pf,{x,v^,t) , 



2 I J ^ 2 

pf,{x,v^,t) p^{x,v^,t) 



(p.) 



d^v^ (fv^ (fv*^ 



XriMvi)) /^,(.^,)ai'^(f,.^,,t)^i^|f^ci\: 



q^,{x,v^,t) 



(Pi) 

y ^^.(^^^^'^(f,^!,^)— ^^^^ 

I 



m,j,q^{x,Vf,,t)d 



(12) 
(13) 

(14) 
(15) 



(16) 



Here, ip^lv^jV^) is some arbitrary function of and -i/^. 

As far as pedestrians of type fi are concerned, we are mostly interested in their density 
(Pn), their mean velocity {vn) and the variance ((5v^,i)^) of their velocity components V/^^i 
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(at a given place x and time t). Since it is formally equivalent and more comparable to 
fluid dynamics, we shall instead search equations for the mass density 

the mean momentum density 

and the mean energy density (in direction i) 



m 



By multiplication of (||) with i'^iv^) = m^, m^v^j, or m^vf^ j2 and integration over 
one obtains the following equations (keeping in mind that the Gaussian surface integrals 
vanish) (see chapter 2.10 of Ref. 



dt 



d 



+ E 



m,. 



((f?M)(^M,a))+QM(l) 



(17a) 
(17b) 



for the mass density, 

^((gM)(^M./3)) 

dt 



777 



(18a) 
(18b) 
(18c) 
(18d) 



for the momentum density, and 



dt 



d 



dx 



fj,,a 



m, 



'''iJ- 

+ E(^m)(^^):^ 

V 

+ E 



m / V 



2 

2 



* 2^ 
fj.,i 



X^iv 



Ml* 



{q,)x7 



(19b) 
(19c) 

(19d) 

(19e) 



for the energy density. Here, we have used the EiNSTElNian summation convention to sum 
over terms in which the Greek indices a, /? or 7 occur twice. 
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3.1 Interpretation 



p,7a| ), ( |18a| ) and ( |19a| ) are the well known hydrodynamic equations (see chapters 2.4 and 



2.10 of Ref. [jl8l). ( |18c|) , ( [L9d|) describe the effects of interactions between two individuals 



of type n and v (for details see section |^). These terms do not vanish as they would do, if 
a conservation of momentum and energy would be fulfilled in a strict sense (see chapter 
2.10 of Ref. |T^|). However, since the individuals try to approach the intended velocity -u^ 



there is a tendency to restore momentum and energy that is described by ( [18b|) , ( |19b|) , 

(0:9^). 



( |17bD , ( |18d| ), ( |19e| ) are additional terms due to individuals who change their type of 
motion. In the following we will assume the special case that these terms as well as the 
terms Qfj,{ipfj,{Vfj,)/m^) due to individuals entering or leaving the system vanish (by 
compensation of inflow into fi and outflow from fi). For concrete situations the quantities 
X'^ and Qfi{.) have to be obtained empirically. 

Pfj,,ai3 is, in thermodynamics, termed the tensor of pressure. Pfj,,ai3n/3 has the meaning of the 
force which is used by the individuals of type n to change their movement when crossing 
a line of unit length / (or, more exactly, the meaning of the component of this force in the 
direction n perpendicular to the line), j^^i is, in thermodynamics, called the heat How. For 
pedestrians it describes the tendency of the velocity variance {{Sv^^iY) to equalize with 
time (see (0)). The variance 

is the thermodynamic equivalent to the absolute temperature T^^i in direction i. Approx- 
imate expressions for p^,q/3 and j^^j are derived in sections 4 and 5. 



3.2 Problems of small densities 

For pedestrian crowds the densities are usually very small. As a consequence, equation 
(§) won't be fulfilled very well and a discrete formulation would be more appealing (see 
||n[). However, we could equivalently start from the continuity equation (PP]), which holds 
better, because the densities are only moderately small. The macroscopic equations 
will be even better fulfilled, because they are only equations for the mean values (^^), 
{v^), (e^,j) and could be also set up by plausibility considerations. 

In order to have small fluctuations of the variables {q^), {v^), (e^^j) with time, in 
equation (|l|) has to be averaged over a finite area A and a finite volume V, which should 
be taken sufficiently large. If T denotes (appart from fluctuations) the time scale for the 
temporal change of {Ofi), (Vf^) and (e^,j), these variables can be also averaged over time 
intervalls At <^ T: 

t+At/2 

1 



{gf,{x,t)){ijf,{x,t)) := — J {g^{x,t')){ij^{x,t')) dt' . 

t~At/2 

Then, (p!7| ) to (p!9| ) will be proper approximations for the movement of pedestrians. 

Another complication by low densities is, that Knudsen corrections have to be taken into 
account [l^: According to these corrections, the "temperature" 9^^i and the tangential 
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velocity {v^^\\) change discontinuously at a boundary dAi, which therefore seems to be 
shifted by a small distance ^ that is comparable to the mean interaction free path (see 
§14 of Ref. ig). 



4 Pedestrians in equilibrium 



In order to calculate p^,Q/3, jfj,,i and x^lu we need the explicit form of (see ([T0|), and 
(|1^) to (p!5D). is the density of individuals of type /i at place x and time t having the 
actual velocity -u^ but arbitrary intended velocity (see (0)). It is directly measurable 
in pedestrian crowds. By integration of (||) over we obtain the theoretical dependence 



dt 



(20) 



J2 C*/^^ + 



(compare to |2^, ^). So, the temporal development of the density is given by a ten- 
dency to walk with the intended velocity (see (||)), by the effects S'^,^ of pair interactions 
(see (^), by the effects C^^, of pedestrians changing their type of motion (see (||)) and by 
the effect of pedestrians entering or leaving the system Ai (see (|ll|), (|l6l)). The last 
two effects shall be be neglected in the following (see the comment in section |3.1|). 



Equation ( pO]) can be solved in a suitable approximation by the recursive method of 
Chapman and Enskog [|, ||. The lowest order approximation presupposes the condi- 
tion dp'^^/dt = of local equilibrium, which is approximately fulfilled by the GAUSSian 
distribution 



1 



-[K,ii-K.,ii))'/(2^?M,ii)+K.±-K,x»V(2e„.x)] 



[21) 



according to empirical data [0, |T^, 



9„ 



< 1 



describes the fact that the velocity variance 6'^,j perpendicular (± 
direction of movement (■u^) is normally less than parallel (||) to it 

For each type p of motion let us perform a particular transformation 



to the mean intended 
1. 



X 



X„ 



X 



X,, 



fu 



8 



Axi] 



P, 



=: e 



/1, 0/3 J 11,13,1 



J, 



fi,a,i 



with 



~p.,al3j 



1 

hi 



This transformation stretches the direction perpendicular to ({/JJ) by the factor 1/6^ and 
simphfies equations ([T7|) to (|l^) to isotropic ones (that means to equations with local 
rotational symmetry). With 



0^, := 61^,1 



we get 

for the pressure and 

for the "heat flow" (see chapter 2.10 of Ref. In addition, the Euler equations 



" fj.,a,i 







dt 

d{V,,p) 
dt 

de^ 

dt 



dt 
dt 



li,a 



1 ^P^,al3 _|_ F^l,|3 



m,. 



dt 



dX„ 



can be derived from (H) (see chapter 16.2 of Ref. pi||). 



(22) 
(23) 

(24) 
(25) 
(26) 



4.1 Behavior on a dance floor 

On a dance floor like that of a discotheque, two types of motion can be found: Type 1 
describing dancing individuals, type 2 describing individuals standing around and look- 
ing on. We can assume to have an isotropic case, that means = 1 and y = 
According to ( |19c| ), the variance {{Sv^ j)^) of the intended velocities f° ^ is causal for the 
temperature 9^^i, that means for the variance ((5f^_j)^) of the actual velocities v^^i. So for 
the temperatures 9i, 6*2 of individuals dancing and individuals standing around it holds 

9i>92, 

since the dancing individuals intend to move with higher variance {{Sv^j)"^) > ((5^2 J^). 
(This is even so in the case of equilibrium, because we have to take the effect of the 
Knudsen corrections into account, see section |3.2| .) The equilibrium condition of equal 
pressure 

pe pe 

-'1 ~ -'2 

now leads to 

6*2 

{Qi) = J^{Q2) < {Q2) 



see (p2|)). Therefore, the dancing individuals will have less density than the individuals 



standing around (see figure |T]). This effect can actually be observed. 
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Figure 1: Behavior on a dance floor: Dancing individuals (filled circles) show a lower 
density than the individuals standing around (empty circles), since they intend to move 
with a greater velocity variance. 

4.2 Propagation of density waves 

In a nearly homogenuous pedestrian crowd with small density variations one can assume 

y^-^xMy^"^^^ WVx.W^O, (l^^)V^/;.^0, F^,;3^0. 
From the EuLER equations ( pi] ) to ( PB]) the equation 



can be derived, then (see chapter 16.2 of Ref. I^Tl)- Subtracting (g^)/6'^x (|26|) from (|2 



and making use of (|22D and (E^), the adiabatic law 



P: 



^ — const. 



can be shown, which leads to the linear wave equation 



^dt^ 

with the adiabatic compressibility 



(see chapter 16.2 of Ref. |^). (For the description of nonlinear waves see |^ and chapter 



2.1 of Ref. [^.) Equation (^) describes the propagation of density waves with velocity 

= (28) 
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On the other hand, the velocity of propagation is given by the mean distance d^ = 1/ ^ (p^) 
of the succeeding individual devided by its mean reaction time C^: 



1 



Inserting this into (pSf), it follows that for small densities (p^) the adiabatic compressibility 
grows with the mean reaction time of individuals according to 



5 Nonequilibrium equations 



In cases of deviations 



5p,iX^, V„ t) := p^(X^, V,, t) - pliX^, V;, t) 



(29) 



from local equilibrium p^(X^, \7^,t) := p^((p/^), (\7^),6'^), we have to find a higher order 
approximation of equation ( pOD than in section ^ If the deviations 5p^ remain small 
compared with p^, we can linearize equation (pO]) around p|^ and get 



E 



(see chapter 15 of Ref. [|T71). Here, r^j, is the mean interaction free time between an 
individual of type p and individuals of type v (see (^5]) and chapter 16.2 of Ref. [pTp. 
From dgl]) and (g2]) to (g6|) one finds [gl|] 



rfp; 



dt 



dpi d{p^) 
d{pf,) dt 

Pi 



+ Pi 



1 



26, 



dt 

- 2 



dOf, dt 



If 



denotes the inverse matrix of 



the relation 



r 



— Sfiv — 



T, 



r, 



dpl 
dt 



E 



r 



tit 



(30) 
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leads, because of PpL = pI^ + Sp^, to the corrected tensor of pressure 

u 

(see flTsD) and the corrected heat flow 



(see (ig)). Here, 



\ dX^^a dX^^a / 



is the tensor of strain, 

are coefficients of the shear viscosity, and 

are coefficients of the thermal conductivity. 

Note, that the effect of restoring the local equilibrium distribution results from the 

tendency to approach the intended velocity distribution p^^{X^, V^, t) with a time constant 
r^, but not from interaction processes as usual (see chapter 13.3 of Ref. ||3^). Therefore 
the viscosity i]^ is dependent on the density (p^) in contrast to ordinary fluids (see pages 
323 and 327 in Ref. [^]). For vanishing densities {p^) — > the interaction rates 1/r^i/ 



become negligible (see (^)) and T^j_y, rj^^, vanish in comparison with r^, i]^, n^, 
respectively. According to (0), the deviation from the local equilibrium distribution 
p^ and, therefore, the viscosity and the thermal conductivity are consequences of finite 



relaxation times r^, r^,y. 



5.1 Effect of viscosity 

For pedestrians the effect of viscosity is not compensated by a gradient of pressure as in 
ordinary fluids, but instead by the tendency of pedestrians to reach their intended velocity 
described by ( |18b|) . In the case of a stationary flow in one direction (that means of one 
type of motion) parallel to the boundaries dM. we have essentially the equation 

= ^MIM = + (d^o^i) - , (32) 

if Tj^ ^ Tj^^ (see (|18a|) , (|18b| )). For a lane of effective width 214^ (with the origin X^^^ = 
in the middle) equation (p2D has the hyperbohc solution 



(^.,11 ) = «||) 
with a boundary layer of width 



cosh(X^,_L/-Dp 



cosh{W/D^ 



(33) 
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In comparison with this, a pressure gradient 

generating the driving force would lead instead to the parabolic solution 

(^M,ll) = ^(W^'-^i±), (34) 

and the mean tangential velocity (V^,||) would depend on the length L of the lane (see 
chapter 3.3 of Ref. [^). The hyperbolic solution ( PB| ) and the parabolic solution ( |51D are 
depicted in figure ||. 



- 



Figure 2: Effect of viscosity (internal friction): Ordinary fluids show a parabolic velocity 
profile (broken line). In contrast to this, a hyperbolic velocity profile is expected for pedes- 
trian crowds (solid line). Whereas in ordinary fluids the internal friction is compensated 
by a pressure gradient, for pedestrian crowds this role is played by the accelerating effect 
of the intended velocity. Due to the Knudsen corrections the fluid slips at the boundary, 
that means the effective width is greater than the actual width. 



6 Effects of interactions 

The scattering rates a*^ of interactions (see (^, (|)) are proportional to the relative 
velocity — "wl/ll of the interacting pedestrians and to the scattering cross section l^y 
(which is a length here and of order of a pedestrian's stride) 
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It is 

1 1 



= {Pu)lA¥,-vA) (35) 

the mean rate of interactions of an individual of type p with individuals of type and 
the corresponding mean interaction free time |jl8|. For the mean relative velocity 



— Vy\\) (see figure H) the following limits can be calculated by making use of (pT]), 
^) and neglecting terms of order 0{5p^): 

I II (^m) - II II II K^l^) - \Vv) II > ^A.> ■ 



Figure 3: The mean relative velocity {\\v^ — Vi,\\) in dependence of \\{v^) — {vy)\\ for the 
special case O^^i = 1 = O^^i. 

Let us introduce 

'^Iv = ^iluii^t^)^ {vu),d,,, du) ■■= r^y{Qu) = T^,ym^{pu) , (37) 



and the total rate of interactions 



^:=E^- (38) 
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Then, 

= e-^*"/^M (39) 

is the probabihty for the possibihty to pass an individual to the right or left, if this would 
need an interaction free time of at least At^ (see chapter 12.1 of Ref. pO|] ). 

k 

is the probability, that two individuals of types \i and v have velocities and Vy before 
their interaction, and the individual of type \x has the velocity v*^ thereafter. We shall 
distinguish three types k of interaction: 

If an individual of type /i is hindered by another individual of type z/, it tries to pass the 
other to the right with probability 'Py.y or to the left with probability 1 — p^j^: 

v*^ = \py,vK^^ - S/s^^v^) + (1 - Pt,u)5iv*^ - S-plv^] 

= S_p^^v^ describes a rotation of velocity to the right side by an angle I3^y in order 
to avoid the hindering pedestrian, v*^ = S^'^^v^ is the inverse rotation to the left side. 

In cases, where it is impossible to avoid the individual of type v having a velocity Vp, the 
individual of type /z tries to walk with velocity v*^ = v,^, if Vi, has a positive component 
■JV ■ > in the intended direction := of motion: 

P^ui^i^, ^v] t^) = (1 - r^)5(t^ - Vp)Q(vp ■ e; > 0) . 

This corresponds to situations, where one moves for a short time within a gap behind a 
pedestrian being in the way (or sometimes, for different directions 7^ el/, in front of it). 
The decision function is defined by 



1 if X is fulfilled 
else. 



If ■ < (that means in the case of a negative component of the hindering pedestrian's 
velocity v,, with respect to the intended direction of movement) it is better for the 
individual (s) to stop (t^ = 0): 



PiM, Vu; t^) = (1 - r,)Si^^ - 0)e(^. ■ e, < 0) . 



This results in 



(n)(n) I ^/^[P/-(^/3..^m) + (1 - Pf.^){Sp^^^f.) - (tv)] , k = l (40.1) 
^^^■{ {l-r,){e,.)[{vy)~{v,)], k = 2 (40.2) 

-(l-r^)(l-(e^,))(t;^). A; = 3 (40.3) 



and 
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— {Qi,){Qu) 



<J* 



k = 1 

k = 2 
k = 3 



-(l-r,)(l-(e,.))«J, 
which allows to calculate ( |18c| ) and ( p.9d|) exphcitly. We have used the abbreviation 



(41.1) 
(41.2) 
(41.3) 



(0M-)(e^, {v,),e,) := {Q{v, ■ > 0)) 



2 



if > 
if y^^ < 



(see figure §) with 

and the GAUSSian error function 



TT 



(1 — r^){l — {Q^u)) is the relative frequency of stopping processes and the relative 
frequency of avoiding processes to the left or right due to interactions. 



6.1 Interpretation 

(a) Development of lanes 

According to (40.1), an asymmetrical avoiding probability p^j, ^ 1 — p^u (see 0) 
leads to a momentum density that tends to the right (for p^j, > 1/2) or to the left 
(for p^i, < 1/2). This momentum density vanishes when the products {q^){Qv) of the 
densities {q^), (Qu) have become zero. Therefore, it causes a separation of different 
types /i 7^ z/ of motion into several lanes (see figure This effect can be observed 
at least for high densities (g^), {Qu) @7 ^ ^ and has the advantage as well as 



the purpose to reduce the total rate 1/f^ of interactions. 

The width of the lanes of two opposing directions 1 and 2 can be calculated from 
the equilibrium condition of equal pressure: 

Pi = P2 , that means (^ii)^i = {§2)02- 

Since in the case of a lane of width and length L consisting of A*"^ individuals 
the relation ^ 

holds for the mass density {g^), we get 

iVi _ fii 
A^2 B2 
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Fi gure 4: The function {Qfxi/} in dependence of 



e^l for the special case 



(i) 



O 



O 



9^0 o • ^ 



o 



o o 



(ii) 




O O i 

Oqo 

Figure 5: (i) Opposite directions of motion normally use separate lanes. Avoiding 
maneuvers are indicated by arrows, (ii) For pedestrians with an opposite direction 
of motion it is advantageous, if both prefer either the right hand side or the left 
hand side when trying to pass each other. Otherwise, they would have to stop in 
order to avoid a collision. The probability p^,, for choosing the right hand side is 
usually greater than the probability (1 —Pnv) for choosing the left hand side. 
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if rriiOi ^ m2^2- That means, the lane width will be proportional to the number 
A*"^ of individuals (see p3[, taking into account the Knudsen corrections described 



in section 3.2 



(b) Crossings 

If the direction fi of motion is crossed by the direction u of motion, it suffers a 
momentum density (40.2), which causes the individuals of type fi to be "pushed" 
partly in direction (tv) of type u. (For the delay effect of crossings see [^.) 

(c) Pedestrian jams 

In order to find out the consequences of (40.3) (that means of stopping processes), 
we can consider the equation 

= -{g,)i{ifi) - {v,)) - (42) 

describing the tendency to walk with the intended velocity t7|J as well as stopping 
processes (see ( |18b|) , (|18cj) ). 



has been introduced for shortness. The stationary solution of equation (|^) is given 

by 



[compare to ||25|). According to ( P5| ) to ( |5D| ) we find 



> 



and, using the abbreviation {bv^^) := — Vy\) , 



^>0, 4((^„,„)=0)=0. (43) 

Due to (^31) a development of pedestrian jams (that means fc^ < 1) is caused by the 
variation 5v^v := — of the velocities. This is even the case for a lane consisting 
of individuals of one type only (where {qu) = for z/ 7^ /x, see (a)), since sj^^ is 
growing with the velocity variance 



ds 



3 



^>o, slie, = o) = o. 



(d) (41.3), describes a loss of variance (a loss of "temperature") by stopping processes. 
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7 Applications 



7.1 Optimal motion 



From sections 5.1 and ^]1| we can conclude that the motion of pedestrians can be optimized 
by 

• avoiding crossings of different directions n of motion, for example by bridges or 
round-about traffic, 

• separation of opposite directions of movement, for example by different lanes for 
each direction (the right line being preferred |25[), or walking at a narrow 
passage by turns, 

• avoiding great velocity variances 6'^, for example by walking in formation (as done 
by the military) p3| . 



• avoidance of obstacles, narrow passages and great densities. 
These rules are applicable to town- and traffic-planning. 



7.2 Maximal diversity of perceptions 

In a museum or super market the individuals will perceive more details (and probably 
buy more goods) if they walk slowly. So the opposite of the rules in 7.1 could be applied 
for planning museums and markets. 



7.3 Critical situations 

In critical situations pedestrians may panic. If the mean total interaction free time (see 
(^)) is less than the mean reaction time (f^, the danger of falling and getting injured is 
great. 

> Cm 

gives a condition for the critical density (p^)^^** of pedestrians that should not be exceeded 
(see (HD to (|8D). 



8 Conclusions 

Starting from the microscopic view of explicit gaskinetic equations we have derived some 
fluid dynamic equations for the movement of pedestrians. These equations are anisotropic 
(that means without local rotational symmetry). They look similar to the equations for 
ordinary fluids, but they are coupled equations for several fluids, that means for several 
types of motion /i, each consisting of individuals having approximately the same intended 
velocity v^. They also contain some additional terms that are characteristic for pedestrian 
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fluids. These terms are due to the tendency of pedestrians to walk with an intended veloc- 
ity, due to pedestrians changing their type (direction) fi of motion, and due to interactions 
between pedestrians (that means due to avoiding maneuvers). 

For high densities (p^) the interactions between pedestrians are very important. As a 
consequence, a development of pedestrian jams and of separate lanes for different direc- 
tions of motion is expected. Pedestrian jams can be understood as an deceleration effect 
due to avoiding maneuvers, and become the worse the greater the velocity variance is. 
The seperation into several lanes is caused by asymmetrical probabilities for avoiding a 
pedestrian to the right or to the left. This asymmetry effect has the advantage to reduce 
the situations, where hindering avoiding behavior is necessary. 

For pedestrian crowds the mechanism of approaching equilibrium is essentially given by 
the tendency to walk with the intended velocity, not by interaction processes as in ordinary 
fluids. As a consequence, the viscosity 77^ (that means the coefficient of internal friction) 
is growing with the pedestrian density. In addition, we have seen that variations within 
pedestrian density will show wave-like propagation with a velocity which depends on 
the mean reaction time. 

Last but not least quantities like "temperature" and "pressure" play another role as in 
ordinary fluids. It can be shown, that the "temperature" (that means the velocity variance) 
6'^ is produced by the variance of the intended velocities, the individuals want to move 
with. As a consequence, two contacting groups of individuals belonging to different types of 
motion can show different "temperatures" . This is the case, for example, on a dance floor. 
On the other hand, whereas a pressure gradient is compensating the effect of internal 
friction in ordinary fluids, for pedestrian crowds this role is played by the accelerating 
effect of the intended velocity. Therefore a hyperbolic stationary velocity profile is found 
instead of a parabolic one. 

9 Outlook 

Current investigations of pedestrian movement deal with the problem how to specify the 
forces fn, and the rates of pedestrians changing their type of motion. These questions 
call for a detailled model for the intentions of pedestrians. 

Pedestrian intentions can be modelled by stochastic laws. They are functions 

• of a pedestrian's demand for certain commodities, 

• of the city center entry points (parking lots, metro stations, bus stops, and so forth), 

• of the location of stores offering the required commodieites, 

• of the expenditures (for example prices, ways), and 

• of unexpected attractions (shop windows, entertainment, and so forth). 

Models of this kind have been developed and empirically tested by BoRGERS and TlM- 
MERMANS 10, A model, which takes into account pedestrian intentions as well as 
gaskinetic aspects will be presented in a forthcoming paper [|lT|. It can be formulated in a 
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way that is also suitable for Monte Carlo simulations of pedestrian dynamics with a com- 
puter. Computer simulations of this kind are an ideal tool for town- and traffic-planning. 
Their results can be directly compared with films of pedestrian crowds (see |TT|). 
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